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Abstract. The dynamics of two interacting spins coupled to separate bosonic baths 
is studied. An analytical solution in Born approximation for arbitrary spectral density 
functions of the bosonic environments is found. It is shown that in the non-Markovian 
cases concurrence "lives" longer or reaches greater values. 
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1. Introduction 

The implementation of more and more efficient nanodevices exploitable in applicative 
contexts like, quantum computers, often requires a highly challenging miniaturizing 
process aimed at packing a huge number of point-like basic elements, whose dynamics 
mimics indeed that of a qubit. Stimulated by such a requirement, over the last few years 
theoretical schemes have been investigated in the language of spin ^ models [1]. Besides 
the simple dynamical behaviour of each elementary constituent these Hamiltonian 
models do indeed capture basic ingredients of several physical situations. In addition, 
spin models allow for the description of the effective interactions in a variety of different 
physical contexts ranging from high energy to nuclear physics |2],[3]. In condensed matter 
physics they capture several aspects of high-temperature superconductors, quantum Hall 
systems, and heavy fermions [H El E] . We point out that Hamiltonians for interacting 
spins can be realized artificially in Josephson junctions arrays [7J or with neutral atoms 
loaded in optical lattices [HI El ED] or still with electrons in quantum dots [H] . 

In this context, a subject deserving a remarkable interest is the entanglement 
dynamics. In view of possible applications it is important to understand the extent 
at which quantum coherences may be protected against the unavoidable degradation of 
the purity of the state in particular when in presence of many-body interactions. 

In this paper we focus our attention on a spin model recently introduced by Quiroga 
[T2] and successively analyzed by other authors [I3]-[T6]. It consists of two interacting 
spins I each one coupled to a separate bosonic bath [12], [13]. Our aim is to study the 
entanglement dynamics of the two spins in the non-Markovian regime. The paper is 
structured as follows. In Section 2 we describe in detail the model. In Section 3 we 
present the analytical solution of the non-Markovian master equation for the reduced 
system constituted by the two interacting spins in the zero temperature limit. In Section 
4 we analyze the entanglement dynamics of the two spins assuming for the environment 
a Lorentz spectral density and an Ohmic spectral density with a Lorentz-Drude cut-off 
function. Finally, conclusive remarks are given in Section 5. 

2. The model 

Our analysis is focused on a system consisting of a pair of interacting spins each one 
coupled to a bosonic bath. The Hamiltonian model that describes such a situation can 
be written as 



contains the free Hamiltonian of the two spins characterized by the same energy e and the 
free Hamiltonian of the two independent bosonic baths characterized by the annihilation 



H = Ho + Hj, 



(1) 



where the unperturbed part Hq 




(2) 



n m 
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and creation operators bni and bl^ [i = 1,2), respectively. As usual af and af are the 
Pauli operators describing the z— th spin {i = 1,2). The interaction Hamiltonian Hj is 
given by 

Hi = K {a^ (T^ + a^a^ ) + ^i^ Qn^iK,! + (^i QnA,! (3) 

n n 
m m 

where K is the strength of the spin-spin interaction and gnA and (7^,2 denote the coupling 
between the spin and its corresponding bosonic bath. In this paper units are chosen 
such that ks = h = 1. In what follows we work in the interaction picture defined by 
Hq, in which the interaction Hamiltonian reads 

Hj{t) = e^^«*i/,e-^^«* = H^'^ + Hf'\t) + i/f ^(t). 



where 



and 



■i(e-uj„,j)t 



(4) 
(5) 
(6) 



In the second order approximation the master equation for the system composed by the 
pair of coupled spins, described by the density operator ps, has the following form [17] 

j^Ps{t) = -z[Hi'\ps{0)]- iyTTTB{[Hj{t),[Hj{T),ps{t)^PB]]). (7) 
Here, we have supposed that p{t) = ps{t) Pb and 



Tr 



B 



Hf'\t),ps{Q)®PB 



0, 



so that 

-li:iB{[Hl{t),ps{^)®PB\)=-l 

Starting from Eq. ([7]) it is possible to obtain 



h\'\ps{^) 



(8) 
(9) 



d 
It 



Ps{t) 



Hf\ps{0)\ + &^ps{t) + £(^^V5(t) + C^'^'^Psit), (10) 



where C^^"* and are the Liouville superoperators associated respectively to the 

interaction Hamiltonian ^ and (jH]) and they are defined by 



H 



H\'\ps{t) 



and 



C^^^^ ps{t) = B'^^\t) [aJps{t),a^]+B^^\t) [oJ.Ps{t)a^] 

+ m{t) [ajps{t),aj]+A^^\t) [aj,ps{t)aj] . 



:iii 



(12) 
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The quantities A^^^{t) and B^^\t) appearing in the previous expression are the so called 
correlation functions, whose explicit form is given by 

A^'\t) = Trfr V|^7.,,f (6l/„,,)i,,e^(^"--)(*--) (13) 

' / . \9n,j\^{K.ibn,j)Bj- 



l 



n '-' 



'^\9n,j\'^{bn,jblj)Bj- 



where {0)Bj = trB^{OpBj}- We point out that the Master equation (fTOj) has been 
obtained under the Born approximation, corresponding to a weak coupling between each 
spin and the relative bosonic bath, truncating the derivation of the master equation to 
the second order of perturbation theory. 



3. Exact solution of the Master equation 



In order to solve the Master equation 
solution for the homogeneous part 
d 



TOl) . it is convenient first step to find a 



(15) 



It is possible and useful for our purpose to separate the equations of motion for the 
diagonal elements of the density operator ps(t) from the ones relative to the off-diagonal 
elements. We have indeed proved that the diagonal elements of psit) have to satisfy the 
following system 

/ Pfit) \ 

Pf{t) 
Pfit) 

\ pfit) I 



d 
dt 



I pm \ 

pfit) 
pfit) 

V pfit) I 



(16) 



where 



( 



\ 



-/3i - /?2 

Pi 





^2 

-2KH - a2 
2KH 



2KH 
-2KH -ai-p2 

^2 



\ 

-ai - a2 j 



(17) 



and 



+ ^. = 5{i)(t) + (18) 

A^^^(t) and B^^\t) being the complex conjugate of A^^\t) and B^^\t), respectively. In 
what follows we will consider the case in which the two bosonic baths are both prepared 
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in a thermal state with T = 0. This assumption in turn imphes that the correlation 
functions reduce to 



A(-^\t)=0, B^^\t) = B{t) = tJ2\9. 



(19) 



Under these hypotheses it is possible to express the solution of the previous system of 
differential equations in the following convenient way 



( Pfit) \ 
Pfit) 
Pfit) 



U^-^\t) 



( p'm \ 

pfio) 

pfiO) 

V P5 (0) J 



(20) 



where 



U^^\t) 



-2G{t) 



-GW(1 



\ 



n _ e-GW)2 



-G(t) 
2 



l-e 



= -G(i) 




(21) 



1 - e-«W 

The function G{t) appearing in the expression for the matrix U'^^\t) is defined as 

G{t) = m + m. 

with 

, 1 - e^(^-^")* + i{e- Un) t 



dTB{T) = ^ \gr, 



(22) 



(23) 



Following an analogous procedure one can also find the off-diagonal matrix elements of 
the density operator psif)- Until now we have obtained a solution for the homogeneous 
part given by Eq. ffTSj) . To find the complete solution of Eq. ffTOj) for the density matrix 
ps{t) we need also a particular solution. In order to do this it is more convenient to 
represent the density matrix ps{t) as a column vector R{t) as follows 

JD(^\ _ ( „ll „22 ^33 44 23 ^32 ^2 34 21 43 ^13 24 31 42 ^14 A1\T /^^x 

- \Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ^ Ps ) y^V 
and to write the solution as 



R(t) = W{t)R{0), 



(25) 



where 



W{0) = he (26) 

is the identity operator of dimension 16 x 16. In Eq. (!25!) W{t) is the propagator that 
has to obey the following equation of motion 
d 



dt 



W{t) = Mo + C'^^'^Wit), 



(27) 
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where Mq is a constant matrix corresponding to the nonhomogeneous term of Eq. fllUI) 
and C^^^ = C^^^\ It is possible to convince oneself that the solution of Eq. fj27|) has the 
following form 

W{t) = V{t) + V{t) [ dTV-\T)Mo, (28) 
Jo 

where V{t) is the propagator associated to the homogeneous equation 

jV{t) = C^^^V{t). (29) 

It is possible to prove that the propagator V{t) has a block diagonal structure 

Vit) = diag (l7W(t),f/f)(t),[/f (t),f/f (t),[/f (t),f/f (t),e-2*W,e-2*W) .(30) 

In the previous expression U^^\t) is the 4 X 4 propagation matrix associated to the 
diagonal elements of ps{t), whereas Uj'^\t) are the 2x2 propagation matrices associated 
to the off-diagonal elements of the ps{t) and they are respectively given by 



U?\t) = '-^{ 1 + 1,,.,. 1 (31) 



and 



u?\t) = — , — I : ). (32) 



At this point we are in condition to explicitly write the density matrix of the two coupled 
spins at a generic time t starting from an arbitrary initial condition. For simplicity, we 
report such a solution in the Appendix. In what follows, instead, we focus on the case in 
which the initial state ofthe pair of coupled spins is the Bell state = ^(|10) — 101)). 
Exploiting the results presented in the Appendix it is possible to demonstrate that the 
state of the reduced system at a generic time instant t can be written in the simple form 

ps{t) = e-^W|^_)(^_| + (1 - e-^W) |00)(00|. (33) 



4. Entanglement dynamics 

As before underlined, the solution we have found has been obtained without specifying 
the spectral properties of the bath. The density matrix ps{t) describing the pair of the 
coupled spins, however, depends on the bath spectral density through the function G{t). 

In this section, exploiting our results, we will analyze some dynamical properties 
of the central system in correspondence to different spectral distributions of the 
environment. In particular, we will examine how the entanglement evolution is affected 
by the choice of the reservoir spectral density. Let's start by considering as first case 
the Lorentzian distribution 
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where 7 and 70 are the reservoir and the system decay rate respectively. This choice in 
turn imphes that the correlation function B{t), as given in the previous section, is 

5(t) = |(l-e-^*) (35) 

and consequently 

G(t) = 7ot+-(e-^*-l). (36) 

7 

We have already demonstrated that starting from the Bell state 

IV^-) = ^(|10)-|01)) (37) 

at a generic time t the density operator describing our system can be written as given 
in Eg. (1551) . It is interesting to analyze how the interaction of the two coupled spins 
with the environments modify the entanglement initially present in the system. To this 
end we consider the time behavior of the concurrence function [TH] of the two spins. 
Exploiting Eq. fl33l) it is easy to demonstrate that in correspondence to any environment 
spectral density, the concurrence function is given by 

C(t) = e-^(*\ (38) 

Thus, when J{uj) assumes the form (IMl) we have 

C{t) = exp f-(7ot + ^(e-^* - 1))^ . (39) 



7 

In Figure [1] we plot C{t) against 7ot for different values of the ratio 7/70. As expected 
in presence of the two baths the concurrence function, starting from its maximum value, 
decreases when the time goes on. However, in the non-Markovian regime, corresponding 
to 7/70 < 2, the entanglement in the two spins remains for longer time with respect to 
the Markovian case. 

Suppose now that the two environments are characterized by an Ohmic spectral 
density with a Lorentz-Drude cut-off function [191 120] 



J{.) = (40) 

where u is the frequency of the bath and Uc is the cut-off frequency. Under this 
hypothesis, putting uq = e/2, the correlation function becomes 

2(^2 

B{t) = -I ^(1 - e-(---^-«)*) (41) 

G{t) = + 4 (.;g - ul)e-^^' sin{u,t) + 



and thus 



+ 8 , ff\,, e~-'' cosM - 8 , ff\,, . (42) 



Exploiting Eq. fl38l) it is possible to analyze the evolution of the degree of entanglement 
of the two spins starting from the Bell state 0371) . The results we have obtained are 




Figure 1. Concurrence function for a Lorentz bath distribution for different values 
of the ratio 7/70 (7/70 — 0-1 (sohd hne), 7/70 = 1 (dashed hne), 7/70 — 10 (dotted 
hne)). The initial state is the Bell state ((37|l . 




Figure 2. Concurrence function for a Lorentz-Drude bath distribution for different 
values of the ratio Wc/wq (wc/^o — 0.1 (solid line), Wc/wq = 1 (dashed line), lOc/ujq = 10 
(dotted line)). The initial state is the Bell state ([57)1 . 
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reported in Figure [2] for different values of the ratio Uc/ujq. Comparing the three plots, 
we may observe that when the spectrum of the reservoir does not completely overlap 
with the frequency of the system, that means Uc ujq, the concurrence function goes 
to zero slower than in the opposite case, ujc ^ ujq. The results we have obtained, 
reported in Figures [U and Ej allow us to note that when the baths are characterized by 
Ohmic spectral densities with a Lorentz-Drude cut-off function, as given in Eq. fHOl) . the 
entanglement initially present in the two spins can be preserved for a longer time with 
respect to the case of a Lorentz bath, at least for some values of the ratio uJc/ujq. 

Following the analysis developed in this section it is also interesting to examine 
the behavior of the system starting from a factorized initial condition instead of an 
entangled one. In what follows in particular we suppose that at t = the two spins are 
in the separable state |1, 0) and we study the time behaviour of the concurrence function 
assuming that the two baths are characterized by the Lorentzian spectral density 
The interaction between the two spins, as expressed by the effective Hamiltonian 
([5]), makes possible the generation of entanglement starting from the factorized initial 
condition before given. On the other hand, in view of the fact that the two spins are 
coupled to two different baths, the quantum correlations that get established in the 
couple of spins will be destroyed. In the non-Markovian regime, however, we expect 
that the entanglement will be preserved for a longer time with respect to the Markovian 
one. This behaviour is confirmed by Figure [3] where we report the time behaviour of the 
concurrence function of the two spins. Looking at this figure we can also observe that 
the degree of entanglement that we can realize in the system starting from the state 
1 1,0) depends on the ratio 7/70- In particular, the maximum value of the concurrence 
function is reached in the non-Markovian case that is 7/70 = 0.1. 

Before concluding we wish to compare our results with the one obtained in the 
Markovian [131 [IB] and post-Markovian [211 122] regimes relatively to the same physical 
system. In order to do this, we concentrate our attention on the temporal behavior of 
the probability Poi(t) of finding the qubit pair in the state |0, 1) supposing that at t = 
the system is prepared in the state |1, 0). In FigureHl where we report Poi(^) in the three 
different regimes, time is scaled in units of the strength K of the spin-spin interaction. 
As shown, when we are in the non-Markovian regime, Poi(^) reaches a maximum value 
that is greater than the one characterizing the Markovian and post-Markovian cases. 
Moreover, as foreseeable in view of the presence of the two baths, in all the regimes the 
probability Poi{t) decays toward zero after reaching its maximum value. However the 
rate of such a decay is smaller in the non-Markovian case. 

5. Conclusions 

In this paper we have analyzed the non-Markovian dynamics of a pair of interacting 
spins coupled to two separate bosonic baths. After deriving the second order master 
equation under the Born approximation only, we have exactly solved it supposing that 
the two bosonic environments are both prepared in a thermal state with T = 0. It 
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Figure 3. Conciirrcncc function for a Lorcntz bath distribution for different values 
of the ratio 7/70 (7/70 = 0.1 (solid line), 7/70 = 1 (dashed line), 7/70 = 10 (dotted 
line)). The initial state is |1,0). 
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Figure 4. Dynamics of the probability to find the system in the state |0, 1) 
(Markovian regime (solid line), post-Markovian regime (dashed line), non-Markovian 
regime (dotted line) for a Lorentz bath distribution with 7/70 = 2). The initial state 
is the separable state |1,0). 
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is important to underline that our solution is valid whatever the initial condition of 
the system as well as the spectral properties of the two baths are. Starting from the 
knowledge of the master equation solution we have studied the temporal behaviour of 
the entanglement established in the couple of interacting spins for different spectral 
densities. The results we have obtained show that in the non-Markovian case the 
concurrence function, adopted in the context of our paper to quantify the degree of 
entanglement in the two spins "lives" or reaches greater values with respect to the 
Markovian regime. Our results motivate further studies on stronger coupling constants 
and non zero temperatures. 
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Appendix 

The full solution for the density matrix of the pair of spins for arbitrary initial conditions 
reads: 



(0), 



p-(t)=e-^W(l 



— e 



«W)p-(0) + 




+ e 



)pf(0) 




2XIm (pf(0)) 



I dre 



2if2r2+G(T) 



pr(t)=e-«W(l-e-«W)p-(0) + 




-2KH'^ 




+2XIm (p|3(0)) 




G(t) 



/ dre 



^r2+G(T) 



pf{t)-i-p'i{t)-pf{t)-pf{t), 
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Jo 

pf{t) = e-«W (Re (pf (0)) + ie-^^'^'lm (pf (0))) 

(pf (0) - p|3(0)) e-^^'^'-G'W /"*dre2^'-'+«(-), 

Jo 

pr(t) = e-^-^W (p^3(0) - e-«Wp^3(o) + ^|4(o)) 

io 

-zii'p^^^Oje-'^-^W-^W Tdre^+^W+^W, 

pr(t) = e-^-*W (p-(0) - e-«Wp-(0) + pr(0)) 

(p^3(0) - pf{0)) e-^^-^W /* cire'^+^W 

Jo 

+zi^p^3^0)e-'^-«W-*W r (ire'^+^M+^W. 
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